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ABSTRACT 

We explore the importance of magnetic-field-oriented thermal conduction in the interaction of su- 
pernova remnant (SNR) shocks with radiative gas clouds and in determining the mass and energy 
exchange between the clouds and the hot surrounding medium. We perform 2.5D MHD simulations 
of a shock impacting on an isolated gas cloud, including anisotropic thermal conduction and radiative 
cooling; we consider the representative case of a Mach 50 shock impacting on a cloud ten-fold denser 
than the ambient medium. We consider different configurations of the ambient magnetic field and 
compare MHD models with or without the thermal conduction. The efficiency of the thermal conduc- 
tion in the presence of magnetic field is, in general, reduced with respect to the unmagnetized case. 
The reduction factor strongly depends on the initial magnetic field orientation, and it is minimum 
when the magnetic field is initially aligned with the direction of shock propagation. The thermal 
conduction contributes to suppress hydrodynamic instabilities, reducing the mass mixing of the cloud 
and preserving the cloud from complete fragmentation. Depending on the magnetic field orientation, 
the heat conduction may determine a significant energy exchange between the cloud and the hot sur- 
rounding medium which, while remaining always at levels less than those in the unmagnetized case, 
leads to a progressive heating and evaporation of the cloud. This additional heating may contrast the 
radiative cooling of some parts of the cloud, preventing the onset of thermal instabilities. 
Subject headings: conduction — magnetohydrodynamics — shock waves — ISM: clouds — ISM: 
magnetic fields — ISM: supernova remnants 



1. INTRODUCTION 

The interaction of the shock waves of supernova rem- 
nants (SNRs) with the magnetized and inhomogeneous 
interstellar medium (ISM) is responsible of the great 
morphological complexity of SNRs and certainly plays a 
major role in determining the exchange of mass, momen- 
tum, and energy between diffuse hot plasma and dense 
clouds or clumps. These exchanges may occur through, 
for example, hydrodynamic ablation and thermal con- 
duction and, among other things, lead to the cloud crush- 
ing and to the reduction of the Jeans mass causing star 
formation. 

The propagation of hot SNR shock fronts in the ISM 
and their interaction with local over-dense gas clouds 
have been investigated with detailed hydrodynamic and 
MHD modeling. The most complete review of this 
problem in the unmagnetized, non-conducting, and non- 
radiative limits is provided by Klein et al. (1994). These 
studies have shown that the cloud is disrupted by the ac- 
tion of both Kelvin-Helmholtz (KH) and Rayleigh- Taylor 
(RT) instabilities after several crushing times, with the 
cloud material expanding and diffusing into the ambient 
medium. An ambient magnetic field can both act as a 
confinement mechanism of the plasma and be modified 
by the interstellar fiow and by local field stretching. Also, 
a strong magnetic field is known to limit hydrodynamic 
instabilities developing during the shock-cloud interac- 
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tion by providing an additional tension at the interface 
between the cloud and the surrounding medium (e.g Mac 
Low et al. 1994; Jones et al. 1996). 

The interaction of the shock with a radiative cloud has 
been only recently analyzed in detail (e.g. Mellema et al. 
2002; Fragile et al. 2004). 2D calculations have shown 
that the effect of the radiative cooling is to break up 
the clouds into numerous dense and cold fragments that 
survive for many dynamical timescales. In the case of 
the interaction between magnetized shocks and radiative 
clouds, the magnetic field may enhance the efficiency of 
the radiative cooling, influencing the size and distribu- 
tion of condensed cooled fragments (Fragile et al. 2005). 

The role played by the thermal conduction during the 
shock-cloud interaction has been less studied so far. In a 
previous paper, Orlando et al. (2005) (hereafter Paper I) 
have addressed this point in the unmagnetized limit. In 
particular, we have investigated the effect of thermal con- 
duction and radiative cooling on the cloud evolution and 
on the mass and energy exchange between the cloud and 
the surrounding medium; we have selected and explored 
two different physical regimes chosen so that either one of 
the processes is dominant. In the case dominated by the 
radiative losses, we have found that the shocked cloud 
fragments into cold, dense, and compact filaments sur- 
rounded by a hot corona which is ablated by the thermal 
conduction. Instead, in the case dominated by thermal 
conduction, the shocked cloud evaporates in a few dy- 
namical timescales. In both cases, we have found that 



2 



Orlando et al. 



the thermal conduction is very effective in suppressing 
the hydrodynaniic instabilities that would develop at the 
cloud boundaries, preserving the cloud from complete de- 
struction. Orlando et al. (2006) and Miceli et al. (2006) 
have studied the observable effects of thermal conduction 
on the evolution of the shocked cloud in the X-ray band. 

Here, we extend the previous studies by investigat- 
ing the effect of the thermal conduction in a magne- 
tized medium, unexplored so far. Of special interest to 
us is to investigate the role of anisotropic thermal con- 
duction - funneled by locally organized magnetic fields - 
in the mass and energy exchange between ISM phases. 
In particular, we aim at addressing the following ques- 
tions: How and under which physical conditions does the 
magnetic-field-oriented thermal conduction influence the 
evolution of the shocked cloud? How do the mass mixing 
of the cloud material and the energy exchange between 
the cloud and the surrounding medium depend on the 
orientation and strength of the magnetic field and on the 
efficiency of the thermal conduction? 

To answer these questions, we take as representative 
the model case of a shock with Mach number = 50 
(corresponding to a post-shock temperature T « 4.7 x 
10^ K for an unperturbed medium with T = 10"' K) 
impacting on an isolated cloud ten-fold denser than the 
ambient medium. Paper I has shown that, in this case, 
the thermal conduction dominates the evolution of the 
shocked cloud in the absence of magnetic field. Around 
this basic configuration, we perform a set of MHD simu- 
lations, with different interstellar magnetic field orienta- 
tions, and compare models calculated with thermal con- 
duction turned either "on" or "off" in order to identify 
its effects on the cloud evolution. 

The paper is organized as follows: in Sect. 2 we de- 
scribe the MHD model and the numerical setup; in Sect. 
3 we discuss the results; and finally in Sect. 4 we draw 
our conclusions. 

2. THE MODEL 

We model the impact of a planar supernova shock 
front onto an isolated gas cloud. The shock propagates 
through a magnetized ambient medium and the cloud is 
assumed to be small compared to the curvature radius 
of the shock"^. The fluid is assumed to be fully ionized 
with a ratio of specific heats 7 = 5/3. The model in- 
cludes radiative cooling, thermal conduction (including 
the effects of heat flux saturation) and resistivity effects. 
The shock-cloud interaction is modeled by solving nu- 
merically the time-dependent non-ideal MHD equations 
(written in non-dimensional conservative form): 



dp 
dt 



V ■ (pu) = 



dpu 



+ V ■ (puu - BB) + VP* = , 



(1) 



(2) 



dpE 



+ V • [u{pE + P*) - B(u • B)] = 

V ■ [B X (77V X B)] - V ■ Fc - nenHA(TX?) 



^ In the case of a small eloud, the SNR does not evolve signif- 
icantly during the shock-cloud interaction, and the assumption of 
a planar shock is justified (see also Klein et al. 1994). 
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V • (uB - Bu) = -V X (r?V x B) , (4) 



where 
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^ = ^+2" +2T 



are the total pressure, and the total gas energy (internal 
energy, e, kinetic energy, and magnetic energy) respec- 
tively, t is the time, p = ^mHriH is the mass density, 
/i = 1.26 is the mean atomic mass (assuming cosmic 
abundances), ttih is the mass of the hydrogen atom, nn 
is the hydrogen number density, u is the gas velocity, T 
is the temperature, B is the magnetic field, rj is the resis- 
tivity according to Spitzer (1962), Fc is the conductive 
flux, and A(T) represents the radiative losses per unit 
emission measure (e.g. Raymond & Smith 1977; Mewe 
et al. 1985; Kaastra & Mewe 2000). We use the ideal gas 
law, P = (7 — l)pe. 

In order to track the original cloud material, we use 
a tracer that is passively advected in the same manner 
as the density. We define Cd the mass fraction of the 
cloud inside the computational cell. The cloud material 
is initialized with Cd = 1, while Cd = in the ambient 
medium^. During the shock-cloud evolution, the cloud 
and the ambient medium mix together, leading to regions 
with < Cci < 1. At any time t the density of cloud 
material in a fluid cell is given by pd = pCd- 

The thermal conductivity in an organized magnetic 
field is known to be highly anisotropic and it can be ex- 
traordinarily reduced in the direction transverse to the 
field. The thermal flux, therefore, is locally split into two 
components, along and across the magnetic field lines, 
Fc = Pii i + F± j, where 
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[9spi]|| [qsa.t]\ 

1 1 



(5) 



[9spi]_L [9sat]± 



to allow for a smooth transition between the classical 
and saturated conduction regime. In Eqs. 5, [qspi]\\ and 
[(7spi]_L represent the classical conductive flux along and 
across the magnetic field lines (Spitzer 1962) 

fepilll = « -5.6 X 10-^TV2 [VT]|| 



[aspi]± = -K±[VT]i « -3.3 X 10- 
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[VT]_ 



(6) 

where [VT]|| and [VT]_l are the thermal gradients along 

and across the magnetic field, and kh and kj_ (in units 
of erg s~^ cm~^) are the thermal conduction coeffi- 
cients along and across the magnetic field lines^, respec- 
tively. The saturated flux along and across the magnetic 
field lines, [qsa.t]\\ and [gsat]^, arc (Cowic & McKec 1977) 

^ We checked that the used numerical scheme ensures that always 
< Cci < 1. 

^ For the values of T, nn and B used here, k^^/k± 10^'' at the 
beginning of the shock-cloud interaction. 
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TABLE 1 

Summary of the initial physical parameters 

rnARACTERizixc THE MIID sI^[^LATI()^'s. 





Temperature 


Density 


Velocity 


ISM 


10* K 


0.1 cm-3 


0.0 


Cloud 


103 K 


1.0 


0.0 


Posf-sh ock 7n r/liurn : 


1.7 X 10" K 


0,1 riii"-^ 


130 km 



feat] 



-sign([VT]||) 5(t>pcl 

-sign([VT]x) 5c^pcl 



(7) 



where Cg is the isothermal sound speed, and is a num- 
ber of the order of unity; we set (j) = 0.3 according 
to the values suggested for a fully ionized cosmic gas: 
0.24 < (j) < 0.35 (Giuliani 1984; Borkowski et al. 1989; 
Fadeyev et al. 2002, and references therein). As dis- 
cussed in Paper I, this choice implies that no thermal 
precursor develops during the shock propagation, consis- 
tent with the fact that no precursor is observed in young 
and middle aged SNRs. 

The initial unperturbed ambient medium is magne- 
tized, isothermal (with temperature Tism = 10"^ K, cor- 
responding to an isothermal sound speed Qsni = 11-5 
km s~^), and uniform (with hydrogen number density 
?^ism = 0.1 cm~^; see Table 1). The gas cloud is in pres- 
sure equilibrium with its surrounding and has a circular 
cross-section with radius rd = 1 pc; its radial density 
distribution is given by 



nc\{r) 



ncio 



cosh [a (t/tc 



(8) 



where ndo is the hydrogen number density at the cloud 
center, r is the radial distance from the cloud center and 
(J = 10. The above distribution describes a thin tran- 
sition layer (~ 0.3 rd) around the cloud that smoothly 
brings the cloud density to the value of the surrounding 
medium^. The initial density contrast between the cloud 
center and the ambient medium is x = '^cio/'^ism = 10. 
The cloud temperature is determined by the pressure bal- 
ance across the cloud boundary. 

The SNR shock front propagates with a velocity w = 
M Cism in the ambient medium, where M. is the shock 
Mach number, and Cism is the sound speed in the inter- 
stellar medium; we consider a shock propagating with 
Ai = 50, i.e. a shock velocity w « 570 km s^^ and a 
temperature Tpsh ~ 4.7 x 10^ K. As discussed in Paper 
I, in this case (for a cloud with r^i = 1 pc and x = 10) 
the cloud dynamics would be dominated by thermal con- 
duction in the absence of magnetic field. The post-shock 
conditions of the ambient medium well before the im- 
pact onto the cloud are given by the strong shock limit 
(Zel'dovich & Raizer 1966). 

Starting from this basic configuration, we consider a set 
of simulations with different initial magnetic field orien- 
tations. We adopt a 2.5D Cartesian coordinate system 
{x,y), implying that the simulated clouds are cylinders 

" A finite transition layer, in general, is expected in real interstel- 
lar clouds due, for instance, to thermal conduction (Balbus 1986; 
see also NaJjamura et al. 2006). 



TABLE 2 

Summary of the MHD simulations. In all runs the shock 
Mach number is = 50, the density contrast is x = 10, 

AND THE CLOUD CRUSHING TIME IS Tec Ri 5.4 X 10^ YR. 



Run 


|B| 
HG 


/3o 


Field 
Comp. 


Therm. 
Cond. 


Rad. 
Losses 


Res."^ 


TSJNT 


u 


oo 




no 


no 




MR 


n 

u 


oo 




no 


yes 




J. i\ 


A 

u 


oo 




yes 


no 






n 
U 


oo 




yes 


yes 




IS iN-JjX'i 




A 
1 


D 


no 


no 




NN-By4 


1.31 


4 




no 


no 


132 


NN-Bz4 


1..S1 


4 


Bl 


no 


no 


132 


1 iN-r>x4 


l.ol 


1 




yes 


no 




TN-By4 


1.31 


4 


By 


yes 


no 


132 


TN-Bz4 


1.31 


4 


Bz 


yes 




132 


TR-Bxl 


2.63 


1 


Bx 


yes 


yes 


132 


TR-Byl 


2.63 


1 


By 


yes 


yes 


132 


TR-Bzl 


2.63 


1 


Bz 


yes 


yes 


132 


TR-Bx4 


1.31 


4 


Bx 


yes 


yes 


132 


TR-By4 


1.31 


4 


By 


yes 


yes 


132 


TR-Bz4 


1.31 


4 


B, 


yes 


yes 


132 


TR-BxlOO 


0.26 


100 


Bx 


yes 


yes 


132 


TR-BylOO 


0.26 


100 


By 


yes 


yes 


132 


TR-BzlOO 


0.26 


100 


B, 


yes 


yes 


132 


TR-Bz4-hr 


1.31 


4 


B, 


yes 


yes 


264 


TR-Bz4-hr2 


1.31 


4 


Bz 


yes 


yes 


528 



''Initial number of zones per cloud radius 

extending infinitely along the z axis perpendicular to the 
{x, y) plane. The primary shock propagates along the y 
axis. In this geometry, we consider three different field 
orientations: 1) parallel to the planar shock and per- 
pendicular to the cylindrical cloud. 2) perpendicular to 
both the shock front and the cloud, and 3) parallel to 
both the shock and the cloud. The magnetic field com- 
ponents along the x and the z axis are enhanced by a 
factor (7 -I- l)/(7 — 1) (where 7 is the ratio of specific 
heats) in the post-shock region (in the strong shock limit; 
Zel'dovich & Raizer 1966), whereas the component along 
the y axis is continuous across the shock. We include runs 
in the strong and weak magnetic field limits, consider- 
ing initial field strengths of |B| = 2.63, 1.31, 0.26, pG 
in the unperturbed ambient medium'', corresponding to 
/3o = 1, 4, 100, 00, where /3o = P/{B^/Sn) is the ratio 
of thermal to magnetic pressure in the pre-shoek region. 
This range of /3o includes typical values inferred for the 
diffuse regions of the ISM (e.g. Mac Low & Klessen 2004) 
and for shock-cloud interaction regions in evolved SNR 
shells (c.g Bocchino et al. 2000). There is no magnetic 
field component c;xclusively associated to the cloud. 

We follow the shock-cloud interaction for 3.5 Tcc, where 
Tec ~ X^^^^c\/w is the cloud crushing time, i.e. the char- 
acteristic time of the shock transmission through the 
cloud; for the conditions considered here (x = 10 and 
A4 = 50), Tcc ~ 5.4 X 10^ yr. Each simulation is repeated 
either with or without thermal conduction for each field 
orientation. Table 2 lists the runs and the initial physical 
parameters of the simulations. 

We solve numerically the set of MHD equations us- 

The unmagnetized case (i.e. |B| = 0) described here is analo- 
gous to the one studied in Paper I except for the fact that in the 
present case the cloud is a cylinder rather than a sphere and has 
smooth boundaries. 
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Fig. 1. — Mass density distribution (gm cm ^) in the (x,y) plane, in log scale, in the simulations NN (left half panels) and TR (right 
half panels), sampled at the labeled times in units of Tcc- The contour encloses the cloud material. 



ing FLASH (Fryxell et al. 2000), a multiphysics code in- 
cluding the PARAMESH library (MacNeice et al. 2000) 
for the adaptive mesh refinement. The MHD equations 
are solved using the flash implementation of the HLLE 
scheme (Einfeldt 1988). The code has been extended 
with additional computational modules to handle the 
radiative losses and the anisotropic thermal conduction 
(see Pagano et al. 2007, for the details of the implemen- 
tation). 

The 2.5D Cartesian (x, y) grid extends between —4 and 
4 pc in the x direction and between —1.4 and 6.6 pc in the 
y direction. Initially the cloud is located at [x, y) = (0, 0) 
and the primary shock front propagates in the direction 
of the y axis. At the coarsest resolution, the adaptive 
mesh algorithm used in the flash code uniformly cov- 
ers the 2.5D computational domain with a mesh of 4^ 
blocks, each with 8^ cells. We allow for 5 levels of refine- 
ment, with resolution increasing twice at each refinement 
level. The refinement criterion adopted (Lohner 1987) 
follows the changes of the density and of the tempera- 
ture. This grid configuration yields an effective resolu- 
tion of ~ 7.6 X 10"'^ pc at the finest level, corresponding 
to K, 132 cells per cloud radius. In Sect. 3.5, we discuss 
the effect of spatial resolution on our results, considering 
the additional runs TR-Bz4-hr and TR-Bz4-hr2 which 
use an identical setup to run TR-Bz4, but with higher 
resolution (w 264 and « 528 cells per cloud radius, re- 
spectively; see Table 2). 

We use a constant inflow boundary condition for the 
post-shock gas at the lower boundary, with free outflow 
elsewhere. For runs with zero magnetic fleld (/3o — oo), 
we use reflecting boundary conditions at x = along the 
symmetry axis of the problem and only evolve half of the 
grid. 

3. RESULTS 

3.1. Dynamical evolution 

Figs. 1 and 2 show the evolution of the mass density in 
the {x, y) plane in the simulations with (3^ — 00 (runs NN 
and TR) and with /3o = 4 (runs NN-Bx4, NN-By4, NN- 



Bz4, TR-Bx4, TR-By4, TR-Bz4). The left (right) half 
panels show the result of models without (with) thermal 
conduction and radiative losses. 

From Fig. 1, we note that the thermal conduction 
drives the cloud evolution in the unmagnetized case 
(/3o = oo; run TR): after the initial compression due 
to the primary shock, the cloud expands and gradually 
evaporates due to the heating driven by the thermal con- 
duction in a few dynamical timescales (see right half pan- 
els in Fig. 1). The heat conduction strongly contrasts the 
radiative cooling of some parts of the cloud and no ther- 
mal and hydrodynamic instabilities (visible in run NN; 
see left-panels in Fig. 1) develops during the cloud evo- 
lution, making the cloud more stable and longer-living 
(the mass mixing is strongly reduced; see Paper I for 
more details). 

We now discuss the effect of the magnetic-field-oriented 
(anisotropic) thermal conduction on the shock-cloud col- 
lision when an ambient magnetic field permeates the 
ISM. We first summarize the expected evolution in the 
presence of an ambient magnetic field, according to the 
well-established results of previous models without ther- 
mal conduction. We distinguish between fields perpen- 
dicular to the cylindrical clouds (i.e. with only and 
By components; referred to as "external" fields by Frag- 
ile et al. 2005) and fields parallel to the cylindrical clouds 
(i.e. with only the B^ component; referred to as "inter- 
nal" fields). In the former case, the magnetic field plays 
a dominant role along the cloud surface and in the wake 
of the cloud where it reaches its highest strength (and 
the plasma (3 its lowest values; e.g. Mac Low et al. 1994; 
Jones et al. 1996). In the case of i?x, the magnetic field 
is trapped at the nose of the cloud, leading to a contin- 
uous increase of the magnetic pressure and field tension 
there (see upper panels in Fig. 2); in the case of By, the 
cloud expansion leads to the increase of magnetic pres- 
sure and field tension laterally to the cloud (see middle 
panels in Fig. 2). In the case of B^ (internal field), the 
magnetic field, being parallel to the cylindrical cloud, 
modifies only the total effective pressure of the plasma 
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Fig. 2. — As in Fig. 1 for tlic simulations with /3o = 4 and tlic magnetic field oriented along x (upper panels), y (middle panels), and z 
(lower panels). The figure shows the distribution in models either without (left half panels) or with (right half panels) thermal conduction 
and radiative losses. For runs NN-Bx4, TR-Bx4, NN-By4 and TR-By4, we plot the magnetic field lines; for runs NN-Bz4 and TR-Bz4, we 
include contours of log(_B^/87r). 
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(Jones et al. 1996); in the case of radiating shocks, the 
additional magnetic pressure may play a crucial role in 
the shocked cloud, preventing further compression of the 
cloud material (Fragile et al. 2005). 

3.1.1. External magnetic fields 

In the case of predominantly external magnetic fields, 
Mac Low et al. (1994) and Jones et al. (1996) have shown 
that the hydrodynamic instabilities can be suppressed 
even in models neglecting the thermal conduction due 
to the tension of the magnetic field lines which maintain 
a more laminar flow around the cloud surface (see also 
Fragile et al. 2005): for a 7 = 5/3 gas, the KH instabil- 
ities are suppressed if (3 < whereas RT instabili- 
ties are suppressed if /3 < (2/7)(x/A1)^ (sec also Chan- 
drasekhar 1961). However, for the parameters used in 
this paper (x = 10 and M. = 50), the magnetic field can- 
not suppress KH instabilities in any of our runs, whereas 
the RT instabilities are suppressed only in runs that lead 
to locally very strong field (/? < 0.05). This can be seen 
in model NN-Bx4 (upper panels in Fig. 2), presenting a 
large field increase at the cloud boundary, compared to 
model NN (Fig. 1): in the latter case the growth of KH 
and RT instabilities at the cloud boundary is much more 
evident than in NN-Bx4. On the other hand, the hy- 
drodynamic instabilities are suppressed more efficiently 
in models including the thermal conduction (runs TR- 
Bx4 and TR-By4) even in cases with low field increase 
(for instance in our By case) as it is evident in Fig. 2 
by comparing models NN-Bx4 and NN-By4 with models 
TR-Bx4 and TR-By4, respectively. 

The thermal exchanges between the cloud and the sur- 
rounding medium strongly depend on the initial field ori- 
entation. Fig. 3 shows the heat flux and magnetic field 
strength distributions in the (.x, y) plane in runs TR- 
Bx4, TR-By4, and TR-Bz4, at time t = 2 t^c- In our 
Bx case (upper panels in Fig. 2), the magnetic field lines 
gradually envelope the cloud, reducing the heat conduc- 
tion through the cloud surface (see left panels in Fig. 3): 
thermal exchanges between the cloud and the surround- 
ing medium arc channelled through small regions located 
at the side of the cloud. The cloud expansion and evap- 
oration are strongly limited by the confining effect of the 
magnetic field (cf. the unmagnetized case TR in Fig. 1 
with model TR-Bx4 in Fig. 2) that becomes up to 30 
times stronger just outside the cloud than inside it (see, 
also, the lower left panel in Fig. 3). The consequent ther- 
mal insulation induces the radiative cooling and conden- 
sation of the plasma into the cloud during the phase of 
cloud compression (t < Tcc). At the end of this phase, 
the cloud material has temperature T « 10^ K and den- 
sity riH ~ 10 cm~^ where primary and reverse shocks 
transmitted into the cloud are colliding; for these values 
of T and nn, the Field length scale (Begelman & McKee 
1990) derived from the ratio of cooling timescale over 
conduction timescale (see Paper I for details) is 

ZwlO^ — w3.2xl0-^pc. (9) 

n-H 

The radiative cooling dominates over the effects of the 
thermal conduction in cold and dense regions with di- 
mensions larger than I. At variance with our unmagne- 
tized case TR, therefore, thermal instabilities develop in 



run TR-Bx4. One of this cold and dense structures is ev- 
ident in Fig. 2 (upper panels) and is located at the cloud 
boundary near the nose of the cloud (at x ~ 0.4 pc and 
y « 3.0 pc) at t = 3 Tcc- 

In the By case, the initial field direction is mostly main- 
tained in the cloud core during the evolution, allowing 
efficient thermal exchange between the core and the hot 
medium upwind of the cloud (see center panels in Fig. 3): 
the core is gradually heated and evaporates in few dy- 
namical timescales. This is illustrated by run TR-By4 
in Fig. 2. On the other hand, the cloud is thermally in- 
sulated laterally where the magnetic field lines prevent 
thermal exchange between the cloud and the surrounding 
medium. Also, a strong magnetic field component along 
the X axis develops in the wake of the cloud and inhibits 
thermal conduction with the medium downwind of the 
cloud. The thermal insulation at the side of the cloud de- 
termines the growth of thermal instabilities where shocks 
transmitted into the cloud collide (see middle panels in 
Fig. 2). 

In both external field configurations, elongated struc- 
tures of strong field concentration are produced on the 
axis downwind of the cloud due to the focalization of 
the magnetized fluid flows there (see upper and middle 
panels of Fig. 2, and lower panels in Fig. 3). These fil- 
amentary structures, identified as "flux ropes" by Mac 
Low et al. (1994), arc formed by magnetic field lines 
stretched around the cloud shape and do not carry a 
significant amount of cloud material (as shown by the 
tracer Cd) although the plasma there moves with the 
cloud (see also Gregori et al. 2000). 

3.1.2. Internal magnetic fields 

Predominantly internal magnetic fields strongly sup- 
press the heat conduction, providing an efficient ther- 
mal insulation of the cloud material (see right panels in 
Fig. 3). In the realistic configuration of an elongated 
cloud with finite length L along the z axis, some heat 
would be conducted along the magnetic field lines. The 
characteristic timescales for the conduction along mag- 
netic field lines is (see Paper I) 

Teond « 2.6 X IQ-'^^ • (10) 

We estimate that the cloud would thermalize in Tcond > 
3.5 Tcc (i-c. the physical time covered by our simula- 
tions), if the length scale of the cloud along the z axis 
is L > 3 pc. In this case, hydrodynamic instabilities de- 
velop at the cloud boundary, being both the magnetic 
field and the thermal conduction not able to suppress 
them. The growth of these instabilities is clearly seen 
in Fig. 2 (lower panels). The combined effect of hy- 
drodynamic instabilities and shocks transmitted into the 
cloud leads to unstable high-density regions at the cloud 
boundaries that trigger the development of thermal in- 
stabilities there (see lower panels in Fig. 2). However, as 
discussed by Fragile et al. (2005), internal magnetic field 
lines are expected to resist comprcission in the shocked 
cloud, thus reducing the cooling efficiency. In fact, in 
our run TR-Bz4, the cloud material is prevented from 
cooling below T w 10^ K. Since the thermal conduction 
does not play any significant role in the shock-cloud in- 
teraction, our Bj. case leads to results similar to those 
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i = l.Q T, 




Fig. 3. — Heat flux (upper panels) and magnetic field strength (lower panels) distributions in the {x, y) plane in the simulations TR-Bx4 
(left panels), TR-By4 (center), and TR-Bz4 (right), at time t = 2 Tec- The arrows in the upper panels describe the heat flux and scale 
linearly with respect to the reference value shown in the upper right corner of each panel. The scale of the magnetic field strength is linear 
and is given by the bar on the right, in units of 10 fj,G. The red contour encloses the cloud material. 

obtained by Fragile et al. (2005) and we do not discuss 
further this case. 



3.2. Role of thermal conduction 

In this section, we study more quantitatively the effect 
of thermal conduction on the cloud evolution and, in par- 
ticular, on the cloud compression and on the magnetic 
field increase. To this end, we use the tracer defined in 
Sect. 2 to identify zones whose content is the original 
cloud material by more than 90%. Then, we define the 
cross-sectional area of cloud material, Ac\{t), as the to- 
tal area in the {x,y) plane occupied by these zones. We 
define the cloud compression (or expansion) as ^ciMcio, 
where A^io is the initial cross-sectional area. We also de- 
fine an average mass-weighted temperature of the cloud 
and an average magnetic field strength associated to the 
cloud as 



{T)c 



Cci pT da 



A(Cci>0.9) 



A(Coi>o.g 



Cci p da 



Cci B da 



A(Cci>0.9) 



(11) 



(12) 



Cci da 

/A(Cci>0.9) 

where we integrate on zones with Cd > 0.9. Note that 
our choice of considering cells with the value of the pas- 
sive tracer Cd > 0.9 is arbitrary. To determine how sen- 
sitive the results are to this value and, in particular, to 
small changes in it, we derive our results also considering 
the values Cd > 0.85 and Cd > 0.95. In all the cases, we 
find that the results derived with the different thresholds 
show the same trend with differences lower than 10%. 
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Fig. 4. — Evolution of cloud compression (upper panels), of average temperature (middle panels), and of average magnetic field strength 
(lower panels) of the cloud for runs which neglect the thermal conduction and the radiation (dot-dashed lines; left panels), for runs which 
include the thermal conduction but neglect the radiation (solid; left panels), for runs which include the radiation but neglect the thermal 
conduction (dotted; right panels) and for runs which include both physical effects (dashed; right panels). The magnetized cases with /3o = 4 
are marked with red (initial magnetic field along x), green (initial B along y) and blue (initial B along z) lines; the unmagnetized cases are 
marked with black lines. The light yellow regions mark the location of solutions which have thermodynamical characteristics in between 
the cases of maximum efficiency of the thermal conduction (models TN and TR) and the cases without thermal conduction (models NN 
and NR). By comparing the position of the magnetized models curves inside the yellow region, it is possible to quantitatively assess the 
degree of suppression of the effects of the thermal conduction by the magnetic fields. 



Fig. 4 shows the cloud compression, v4ci/^cio, the av- 
erage temperature of the cloud, (T)ci, normalized to 
the post-shock temperature of the surrounding medium 
(Jpsh = 4.7 X 10^ K), and the average magnetic field 
strength associated to the cloud, {B)c\, normalized to 
the initial field strength {Bq = 1.31 /iG, corresponding 
to /3o = 4) as a function of time for models neglecting 
thermal conduction and radiation (hereafter NNs mod- 
els), for models including conduction but neglecting radi- 
ation (TNs models), and for models including both con- 
duction and radiation (TRs models); we also include the 
results derived from the unmagnetized case NR with ra- 
diative cooling and without thermal conduction. The 
figure shows both the magnetized cases with (3q = A and 



the unmagnetized cases (see Table 2). 

In aU the NNs models either with (NN-Bx4, NN-By4, 
and NN-Bz4) or without (NN) the magnetic field, the 
evolution of the cloud compression and of the average 
cloud temperature is roughly the same (see left panels in 
Fig. 4). The cloud is initially compressed over a timescale 
t w Tec due to the ambient post-shock pressure; during 
this phase (T)ci rapidly increases. After t « Tcc, the 
cloud partially reexpands, leading to a decrease of (T)ci. 
In the last phase (i > 2.0 Tcc), the cloud is compressed 
again by the interaction with the "Mach stem" formed 
during the reflection of the primary shock at the sym- 
metry axis, and (r)ci increases; later Ad/ Ado continues 
to decrease, because of the mixing of the cloud material 
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with the ambient medium (see Sect. 3.3; see also Paper 
I), while (T}ci stabilizes at w 0.17 Tpsh- 

The field increase in the cloud material depends on the 
initial configuration of B (see lower left panel in Fig. 4). 
In the case of external fields {B^ and By components), B 
is mainly intensified due to stretching of field lines due 
to sheared motion. In the _Bx case, the magnetic field 
undergoes the greatest increase and (i?)ci keeps increas- 
ing during the whole evolution. In fact the field is mainly 
intensified at the nose of the cloud where the background 
flow continues to stretch the field lines during the evolu- 
tion (see upper panels in Fig. 2). In the By case, the field 
increase occurs mainly at the side of the cloud where the 
field lines are stretched along the cloud surface. In the 
case of internal fields {B^ component), the field increase 
is due to squeezing of field lines through compression. 
{B)ci: therefore, follows the changes in Ad/Ado, since 
the field is locked within the cloud material. Thus the 
greatest field increase occurs at t ~ Tcc when the shocks 
transmitted into the cloud collide. 

The effects of thermal conduction are greatest in the 
unmagnetized model (TN) which can be considered an 
extreme limit case (see left panels in Fig. 4). During 
the first stage of evolution {t < 0.8 Tcc), the cloud is 
heated efficiently by the thermal conduction and its av- 
erage temperature increases rapidly to ~ 0.5 Tpsh- As 
a consequence, the prcssiirc inside the cloud increases 
and the cloud reexpands earlier than in model NN. Af- 
terwards, the average cloud temperature, (T)ci, keeps 
increasing up to ~ 0.9 Tpsh at t = 3.5 Tcc- 

In the case of predominantly external magnetic fields 
(models TN-Bx4 and TN-By4), the thermal conduction 
still plays a significant role in the cloud evolution, al- 
though its effects are not as large as in the mnnagnetized 
case (TN). During the initial compression, the thermal 
conduction contributes to the cloud heating: the aver- 
age temperature of the cloud reaches values larger than 
in models neglecting the conduction (compare TNs with 
NNs models in the left panels of Fig. 4). This effect is 
greatest in the By case which is the configuration of field 
lines that allows the most efficient thermal exchange be- 
tween the cloud and the hot environment (see Sect. 3.1). 
At t = 3.5 Tcc, (T)ci in TNs models reaches values larger 
than in NNs models (w 0.21 Tpsh in the B^ case and 
« 0.33 Tpsh in the By case). For internal magnetic fields, 
the thermal conduction plays no role in the evolution of 
the shocked cloud, being strongly ineffective due to B 
(see Sect. 3.1). As a consequence, the TN-Bz4 model 
leads to the same results as NNs models. 

In general, therefore, the effects of the thermal con- 
duction in the presence of an ambient magnetic field 
are reduced with respect to the corresponding unmag- 
netized case, but not entirely suppressed. This can be 
seen in Fig. 4, where we have marked in light yellow the 
region between the fully conductive unmagnetized case 
(TN) and the case without thermal conduction (NN). 
The magnetized TNs models are always within this re- 
gion, meaning that the effects of the thermal conduction 
are never as large as in the unmagnetized case (TN) but 
not completely suppressed as in the model NN. 

We also note that the thermal conduction influences 
indirectly the magnetic field increase. The main changes 
are in the By case and are due to the larger expansion 
of the cloud that reduces the increase of the field associ- 



ated to the cloud, being the field locked within the cloud 
material. 

In our unmagnetized case TR (including thermal con- 
duction and radiative cooling), the thermal conduction 
prevents the onset of thermal instabilities, and the evo- 
lution of the shocked cloud is the same as found in the 
TN model. At variance with our mimagnetized case TR, 
Fig. 4 shows that thermal instabilities develop in all our 
magnetized TRs runs, being the effects of thermal con- 
duction reduced by the magnetic field- The effects of 
radiative cooling are very strong for internal fields (our 
Bz case; see run TR-Bz4 in Fig. 4). In this case, the 
heat conduction is totally suppressed by the magnetic 
field and the evolution of Ad/Ado and of (T)ci are the 
same as those found in the unmagnetized case with ra- 
diative cooling and without thermal conduction (model 
NR); at t = 3.5 Tcc, nm TR-Bz4 (and NR) shows the 
largest cloud compression {Ad /Ado ~ 0.1) and the low- 
est cloud average temperature ((T)ci ~ 0.12 Tpsh)- In 
the case of external fields (runs TR-Bx4 and TR-By4), 
the effects of heat conduction are reduced but not sup- 
pressed and the results are intermediate between those 
derived for runs NR and TR (i.e. within the light yellow 
region in the right panels in Fig. 4). The cooling effi- 
ciency is largely reduced in our By case (run TR-By4), 
namely that with the magnetic field configuration that 
allows the most effective thermal conduction. 

3.3. Mass mixing and energy exchange 

We use the tracer to derive the cloud mass, Md, as the 
total mass in zones whose content is the original cloud 
material by more than 90%, 

Md = L [ Cd p da , (13) 

Ja{Cci>0.9) 

where L is the cloud length along the z axis, and the 
integral is done on zones with Cd > 0.9. We investigate 
the mixing of cloud material with the ambient medium 
by defining the remaining cloud mass as Md/Mdo, where 
Mcio is the initial cloud mass. 

The tracer allows us to investigate also the energy ex- 
change between the cloud and the surrounding medium; 
we derive the internal energy, Xd, and the kinetic energy, 
/Ccl of the cloud as 

Tci = T / Cdpeda, (14) 

Ja{Cci>0.9) 

^ci = ^ / C'ci p|u|^ da , (15) 

^ JA(Cci>0.9) 

where again L is the cloud length along the z axis, and 
the integral is done on zones with Cd > 0.9. We also 
define the total energy of the cloud as 

Ed = Id + ICd - (16) 

Fig. 5 shows the evolution of the cloud mass, Md/Mdo, 
for NNs, TNs, and TRs models; again we also include 
the unmagnetized case with radiative cooling and with- 
out thermal conduction (model NR). Both unmagnetized 
cases and magnetized cases with /?o = 4 are shown. In 
models without thermal conduction and radiation (NNs 
models), the hydrodynamic instabilities drive the mass 
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Fig. 5. — Presentation as in Fig. 4 for the evolution of the cloud mass (upper panels), of the internal energy of the cloud (middle panels), 
and of the kinetic energy of the cloud (lower panels). 



mixing of the cloud®. The mass loss rate of the cloud, 
rhci, increases significantly after 1.5 Tcc (i.e. after the 
hydrodynamic instabilities have fully developed at the 
cloud boundary), with rhd ~ 1.5 x 10~^Lpc Mq yr~^, 
where Lpc is the cloud length along the z axis in units of 
pc: ~ 20% of the cloud mass is contained in mixed zones 
at t = 3.5 Tcc. The only exception is run NN-Bx4 (~ 15% 
of the cloud mass is in mixed zones at t = 3.5 Tcc), being 
in this case RT instabilities partially suppressed by the 
magnetic field (compare run NN-Bx4 with runs NN-By4 
and NN-Bz4 in Fig. 2). 

In TNs models with external magnetic fields (TN- 
Bx4 and TN-By4), the mass loss rate of the cloud 
is less efficient than in NNs models with md 6 x 
10~^Lpc Mq yr^^ (^ 10% of the cloud mass is in mixed 
zones at t = 3.5 Tcc). In fact, in these cases the thermal 

* This is also true in our magnetized cases because, for the pa- 
rameters used in this paper {A4 = 50 and x = 10) i the hydrody- 
namic instabilities are partially suppressed by the magnetic field 
only in runs evolving to strong fields (see Sect. 3.1). 



conduction suppresses most of the hydrodynamic insta- 
bilities and the mass loss mainly comes from the cloud 
evaporation driven by the thermal conduction rather 
than from hydrodynamic ablation. Note that our unmag- 
netized TN model is an extreme limit case in which the 
hydrodynamic instabilities are totally suppressed by the 
thermal conduction which drives the cloud mixing; in this 
case, the mass loss rate is md ~ 1.5 x 10~^Lpc Mq yr~^ 
(~ 5% of the cloud mass is in mixed zones att = 3.5 Tcc). 

In magnetized TRs models, the onset of thermal in- 
stabilities increases the mass loss rate of the cloud with 
respect to the unmagnetized case (ttIci ranges between 
1.5 X lO-^Lpc Mq yr^i and 4 x lO'^Lpc Mq yr"!) 
due to the fragmentation of the cloud in dense and 
cold cloudlets. We expect, therefore, that the larger 
the amount of cloud mass mixed with the surrounding 
medium at the end of the evolution, the more limited the 
thermal exchange between the cloud and the hot ambi- 
ent medium (and, therefore, the greater the efficiency of 
radiative cooling). In fact, the upper right panel in Fig. 5 
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shows that the mass mixing has the greatest efficiency in 
run TR-Bz4 (i.e. in the case with the thermal conduction 
totally suppressed) which shows a mass loss rate of the 
cloud similar to that derived from the unmagnetized NR 
model. On the other hand, in runs TR-Bx4 and TR-By4, 
the mass mixing is intermediate between those derived 
with runs NR and TR. 

Fig. 5 also shows the evolution of internal (middle 
panels) and kinetic (lower panels) energy of the cloud, 
normalized to the initial total energy of the cloud, Edo- 
Among the magnetized cases considered, the greatest val- 
ues of Tc\ are reached in our By case which is the field 
configuration that allows the most efficient thermal ex- 
change between the cloud and the environment; the in- 
crease of Xci is due to the heat conducted to the shocked 
cloud. Also, the By case leads to the greatest values of 
/Cci because the cloud has a larger cross-sectional area 
(because of the larger cloud expansion due to the heat- 
ing driven by heat conduction; see upper panels in Fig. 4) 
and offers, therefore, a larger surface to the pressure of 
the shock front responsible of the cloud acceleration. 

3.4. Role of the initial field strength 

In this section we explore the effects of the initial field 
strength on the mass mixing and energy exchange of the 
cloud. Fig. 6 shows the evolution of the cloud mass, 
Mci/Mcio (upper panel), and of the total (internal plus 
kinetic) energy of the cloud, Ec\/Ec\o (lower panel), for 
magnetized TRs models with different values of (3q. We 
discuss here only the cases of predominantly external 
magnetic fields {B^ or By case) since no significant de- 
pendence on the initial field strength has been found in 
the case of predominantly internal magnetic fields {B^ 
case). 

Fig. 6 shows that the initial field strength plays a sig- 
nificant role in the B^ case. In particular, models with 
greater values of po show a more efficient mixing of the 
cloud material and a less rapid increase of the cloud en- 
ergy. As discussed in Sect. 3.3, in the B^ case, the rate 
of mass-loss from the cloud is mainly driven by abla- 
tion through the hydrodynamic instabilities (being the 
thermal conduction strongly suppressed by the magnetic 
field). On the other hand, in the case of external fields, 
the instabilities can be dumped by the magnetic field, de- 
pending on its strength (see Sect. 3.1.1). For instance, 
in the i?x case with Pq — 4, we found that the RT insta- 
bilities are mostly suppressed by the magnetic field (see 
upper panels in Fig. 2). In the B^ case with Pq = 100, 
instead the magnetic field is too weak to dump the hy- 
drodynamic instabilities over the timescales considered; 
these instabilities, in turn, lead to the formation of re- 
gions dominated by the radiative cooling, triggering the 
development of thermal instabilities. Both the hydrody- 
namic and the thermal instabilities determine the cloud 
mass mixing (which is higher for higher values of Pq). In 
addition, the thermal instabilities reduce the increase of 
the cloud energy (which is less rapid for higher Pq) due 
to significant radiative losses. 

In the By case, the initial field strength has a smaller 
influence on the dynamic and thermal evolution of the 
cloud than in the B^ case (see Fig. 6). In addition, at 
variance with the B^ case, models with greater values 
of /3o show a less efficient mixing of the cloud material 
and a more rapid increase of the cloud energy. In the By 




Fig. 6. — Evolution of the cloud mass (upper panel) and of 
the total energy of the cloud (internal plus kinetic; lower panel) 
for runs including both the thermal conduction and the radiative 
cooling (TRs models). The figure shows the simulations with the 
magnetic field oriented along x (red lines) or y (green) and with 
/3o = 1 (dotted lines), 4 (solid), and 100 (dashed). 

case, in fact, the hydrodynamic instabilities responsible 
of the mass mixing are mainly suppressed by the thermal 
conduction rather than by the magnetic field as in the i?x 
case. As a consequence, the higher the value of /So, the 
more effective the thermal conduction in suppressing the 
instabilities and in heating the plasma, the less efficient 
the cloud mass mixing and the more rapid the increase 
of the cloud energy. 

3.5. Effect of spatial resolution 

The effective resolution adopted in our simulations is 
~ 132 cells per cloud radius, a value above the resolu- 
tion requirements suggested by Klein et al. (1994) for 
non-radiative clouds. However, for radiative clouds, we 
expect that the details of the plasma radiative cooling 
depend on the numerical resolution: a higher resolution 
may lead to different peak density and hence influence 
the cooling efficiency of the gas, preventing further com- 
pression of the cloud. In the non-conducting regime. 
Fragile et al. (2005) found that the results generally con- 
verge for simulations with resolution larger than 100 cells 
per cloud radius ( ^10% differences). In the simula- 
tions presented here, the thermal conduction partially 
contrasts the radiative cooling in the case of external 
fields (-Bx or -By), alleviating the problem of numerical 
resolution (see also Paper I). 
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Fig. 7. — Presentation as in Fig. 6 for runs TR-Bz4 (solid lines), 
TR-Bz4-lir (dashed), and TR-Bz4-hr2 (dotted). 

In order to check if our adopted resolution is sufficient 
to capture the basic cloud evolution over the time inter- 
val considered, we compare three simulations (TR-Bz4, 
TR-Bz4-hr, and TR-Bz4-hr2) with different spatial reso- 
lution (132, 264, and 528 zones per cloud radius, respec- 
tively) for the Bz case with /3 = 4, namely one of the cases 
in which the growth of hydrodynamic and thermal insta- 
bilities is most prominent and the effect of thermal con- 
duction (contrasting the development of hydrodynamic 
and thermal instabilities) is negligible. Since this case 
is one of the most demanding for resolution, it can be 
considered a worst case comparison of convergence. 

Figure 7 compares the evolution of the cloud mass, 
Mci/Mcio, and of the total energy of the cloud, Ed/Edo, 
for the three simulations TR-Bz4, TR-Bz4-hr, and TR- 
Bz4-hr2. In general, we find that the results obtained 
with the three simulations agree quite well in their quali- 
tative behavior, showing differences < 10%. In runs TR- 
Bz4-hr and TR-Bz4-hr2, the remaining cloud mass and 
the total energy of the cloud are, in general, systemati- 
cally higher than in run TR-Bz4. The larger mass mixing 
in TR-Bz4 is driven by the higher diffusion of the low- 
resolution grid down to the very small structures which 
tend to smear out concentrated density peaks, promot- 
ing mass mixing. The slightly lower energy of the cloud 
in TR-Bz4 is a consequence of the larger mass mixing 
derived in this run with respect to TR-Bz4-hr and TR- 
Bz4-hr2. Note that, in runs showing the onset of ther- 
mal instabilities (i.e NRs and TRs models), the size of 



the latters reaches the resolution limit toward the end of 
the simulations when the relevant physical processes are 
already at a late stage. 

4. SUMMARY AND CONCLUSION 

We investigated the importance of magnetic-field- 
oriented thermal conduction in the interaction between 
an isolated elongated dense cloud and an interstellar 
shock-wave of an evolved SNR shell through numerical 
MHD simulations. To our knowledge, these simulations 
represent the first attempt to model the shock-cloud in- 
teraction that simultaneously considers magnetic fields, 
radiative cooling, and anisotropic thermal conduction. 
Our findings lead to several conclusions: 

1. In general, we found that the effects of thermal con- 
duction on the evolution of the shocked cloud are 
reduced in the presence of an ambient magnetic 
field with respect to the unmagnetized cases in- 
vestigated in Paper I. The efficiency of anisotropic 
thermal conduction strongly depends on the ini- 
tial magnetic field orientation and configuration. 
This efficiency is the largest when the initial B 
is aligned with the direction of propagation of the 
shock front, and is the smallest when B is aligned 
with the cylindrical cloud, namely when the heat 
conduction is completely suppressed by the mag- 
netic field. 

2. We found that the hydrodynamic instabilities are 
suppressed efficiently by the anisotropic thermal 
conduction when the initial magnetic field is per- 
pendicular to the cylindrical cloud (a configuration 
referred to as "external fields" ). On the contrary, in 
the case of B parallel to the cylindrical axis of the 
cloud (i.e. when the field has component only along 
the z axis - internal field), hydrodynamic instabili- 
ties develop at the cloud boundary. We found that, 
for the parameters of the simulations chosen, the 
magnetic tension is unable to suppress alone the 
hydrodynamic instabilities. 

3. As for thermal instabilities, we found that, depend- 
ing on the magnetic field orientation, the heat flux 
contributes to the heating of some parts of the 
cloud, reducing the efficiency of radiative cooling 
there, and preventing any thermal instability. 

4. The mass loss of the cloud due to mixing with 
the surrounding medium is mainly driven by hy- 
drodynamic instabilities; in the case of external 
fields (initial B perpendicular to the cylindrical 
cloud) the anisotropic thermal conduction reduces 
the mass mixing of the cloud. In any case, the mass 
loss rate is larger than that in the corresponding 
unmagnetized case {rhd ~ 1.5 x 10~''Lpc Mq yr~^, 
i.e. ~ 5% of the cloud mass is in mixed zones 
at t = 3.5 Tec), but can get very high when 
the thermal conduction is completely suppressed 
(mci ^ Ax lO-'^Lpc Mq yr-i, i.e. ~ 45% of the 
cloud mass is in mixed zones at t = 3.5 t^c)- 

5. The thermal conduction mostly rules the en- 
ergy exchange between the cloud and surrounding 
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medium. The exchange is favored when the mag- 
netic field configuration is such that the conductive 
flow is not suppressed (i.e. external field configu- 
rations, Bx and By cases), but it is never as high 
as in the absence of magnetic field. In the By case, 
the cloud core is efiRciently heated and evaporates 
in few dynamical timescales. 

6. In general, the initial magnetic field strength has a 
small influence on the dynamic and thermal evolu- 
tion of the shocked cloud for the ranges of values 
explored in this paper (namely 0.26 jiG < |B| < 
2.63 iiG). 

It is worth noting that some details of our simulations 
depend on the choice of the model parameters. For in- 
stance, the onset of thermal instabilities or the evapora- 
tion of the whole cloud depends on the initial shock Mach 
number, and on the density and dimensions of the cloud. 
The cases that we present here (i.e. A4 = 50, x = lOj 
and different configurations of B) are representative of 
a regime in which both the thermal conduction and the 
radiative cooling play an important role in the evolution 
of the shocked cloud. Nevertheless, our analysis proves 
that anisotropic thermal conduction can not be neglected 
in investigations of the evolution of shocked interstellar 
clouds. 

In our simulations, we consider laminar thermal con- 
duction, although regions of strong turbulence of differ- 
ent strength and extent develop in the system (for in- 
stance, at the shear layers along the cloud boundary or 
at the vortex sheets in the cloud wake). In fact, the 
turbulence in these regions may have a significant ef- 
fect on thermal conduction, leading to significant devi- 
ations of thermal conductivity from its laminar values 
(e.g. Narayan & Medvedev 2001; Lazarian 2006); in some 
cases, the turbulence may enhance the heat transfer, ex- 
ceeding the classical Spitzer value (Lazarian 2006). As a 
result, thermal conduction may be not only anisotropic 
(in the presence of the magnetic field) but also "inhomo- 
geneous" due to the presence of turbulence. However, 
even modeling accurately the turbulent thermal conduc- 
tivity, we do not expect significant changes in the results 
of our Bz case, being the thermal conduction strongly 
ineffective in the whole spatial domain; in the remaining 
cases (i?x and By), our modeled thermal conductivity 
could be underestimated in regions of strong turbulence, 
affecting some details of the simulations but not the main 
conclusion of the paper that, in general, anisotropic ther- 
mal conduction can play an important role in the evolu- 
tion of the shocked cloud. 

Note also that the field configurations studied in this 
work are highly idealized. More realistic fields are ex- 
pected to have more complex topologies and, often, the 
field can be tangled and chaotic. In the latter case, the 
thermal conduction will approach isotropy, whereas the 
effect of MHD turbulence is expected to partially sup- 
press the heat transfer within a factor ~ 5 below the 
classical Spitzer estimate^ (Narayan & Medvedev 2001; 
Lazarian 2006). The shock-cloud collision in the presence 
of an organized ambient magnetic field, discussed here, 

^ As already discussed, the MHD turbulence can even enhance 
the heat transfer in some cases (see Lazarian 2006). 



and that in the absence of magnetic field can be con- 
sidered as extreme cases: the former leading to highly 
anisotropic thermal conduction, the latter to the clas- 
sical Spitzer thermal conduction. The case of chaotic 
magnetic field is expected to fall in between these two. 

Our simulations were carried out in 2.5D Cartesian 
geometry, implying that the modeled clouds are elon- 
gated along the z axis. This choice is expected to affect 
some details of the simulations but not our main con- 
clusions. Adopting a 3D Cartesian geometry and model- 
ing a spherical cloud, the highly symmetric shock trans- 
mitted into the cloud converging on the symmetry axis 
would lead to compression stronger than those found in 
our 2.5D simulations, enhancing the radiative cooling. 
Also, 3D simulations would provide an additional degree 
of freedom for hydrodynamic instabilities, increasing the 
mass loss rate of the cloud in the eases in which the mass 
mixing of cloud material is driven by instabilities. Note 
that, for a spherical cloud, our and B^ cases no longer 
differ. 

Finally we assume, in our simulations, that the cloud 
and the ambient material have the same composition, 
implying that microscopic mass mixing due to shear in- 
stabilities would be irrelevant. In a more realistic condi- 
tion, a cold dense cloud may have a different composition 
from the hot ambient flow and the degree of microscopic 
mixing may translate into different spectral signatures of 
the system. In this case, species diffusion could also be 
important, along with thermal conduction, to determine 
the degree of microscopic mixing of the materials and, 
consequently, one would have to ask about the typical 
values of the Lewis number (i.e. the ratio of thermal 
diffusivity to mass diffusivity) in the system. 

It is worth emphasizing that the quantitative results 
of our simulations depend on the physical parameters 
of the model (shock Mach number, density contrast and 
dimension of the cloud, etc.) as well as on the basic 
assumptions of the model (geometry of the cloud, ge- 
ometry of the ambient magnetic fleld, laminar thermal 
conduction, composition of the cloud and of the ambient 
medium, etc.). Nevertheless, our results undoubtedly 
show that the magnetic-field-oriented thermal conduc- 
tion can play an important role in the evolution of the 
shock-cloud interaction (which depends on the magnetic 
field orientation and configuration) and, in particular, in 
the mass and energy exchange between the cloud and 
the hot surrounding medium. We conclude, therefore, 
that a self-consistent and quantitative description of the 
interaction between magnetized shock-waves and inter- 
stellar gas clouds should include the effects of thermal 
conduction. 

The results presented here are interesting for the study 
of middle-aged SNR shells expanding into a magnetized 
ISM and whose morphology is affected by ISM inhomo- 
geneities (for instance, G272.2-3.2, e.g. Egger et al. 1996; 
Cygnus Loop, e.g. Patnaude et al. 2002; Vela SNR, e.g. 
Miceli et al. 2005). It will be further interesting to ex- 
tend the present study, by modeling the shock-cloud in- 
teraction in 3D with radiative cooling, anisotropic ther- 
mal conduction, and magnetic field included and, even, 
considering detailed comparisons of model results with 
observations. 
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